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Abstract. We briefly describe the three existing scenarios for forming 
massive stars and emphasize that the arguments often used to reject the 
accretion scenario for massive stars are misleading. It is usually not 
accounted for the fact that the turbulent pressure associated to large 
turbulent velocities in clouds necessarily imply relatively high accretion 
rates for massive stars. 

We show the basic difference between the formation of low and high 
mass stars based on the values of the free fall time and of the Kelvin- 
Helmoltz timescale, and define the concept of birthline for massive stars. 

Due to D-burning, the radius and location of the birthline in the HR 
diagram, as well as the lifetimes are very sensitive to the accretion rate 
M accr . If a form M accr oc A(M/Mq) </3 is adopted, the observations in 
the HR diagram and the lifetimes support a value of A ~ 10 -5 Mq • yr _1 
and a value of <p ;> 1. Remarkably, such a law is consistent with the 
relation found by Churchwell (1998) and Henning et al. (2000) between 
the outflow rates and the luminosities of ultra-compact HII regions, if we 
assume that a fraction 0.15 to 0.3 of the global inflow is accreted. The 
above relation implies high M accr ~ 10 -3 Mq • yr _1 for the most massive 
stars. The physical possibility of such high M accr is supported by current 
numerical models. 

Finally, we give simple analytical arguments in favour of the growth 
of M accr with the already accreted mass. We also suggest that due to 
Bondi-Hoyle accretion, the formation of binary stars is largely favoured 
among massive stars in the accretion scenario. 



1. Introduction 

Numerous new observations in radio, IR, optical, UV and X-rays, together with 
many numerical models, have contributed to the progress of the field of star 
formation. However, the formation of massive stars is still a major unsolved 
problem in stellar evolution. 

The answer brought to the question "How do massive stars form ?" has a 
major impact not only for stellar evolution, but also for spectral and chemical 
evolution of galaxies as well as for cosmology. Several competing processes have 
now been identified leading to different scenarios for the formation of massive 
stars. Nevertheless, this does not mean that the authors of these scenarios are 
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competing. Rather they all look for the proper answer to the above major 
question. 

2. The various scenarios for the formation of massive stars 

One can identify three different scenarios: 

a) The classical scenario. 

This is the pre-MS evolution at constant mass with bluewards horizontal tracks 
in the HR diagram, moving from the Hayashi line to the zero age main sequence 
(ZAMS), (c/. Iben 1963). The timescale is the Kelvin-Helmholtz timescale 
^kh ~ G M 2 /RL which is about 0.5 to 2% of the MS lifetime. For example, £kh 
is about 3 TO 4 yr for a 30 Mq star. Due to lots of evidence of mass accretion (see 
below, § |3]), this scenario is no longer supported, although it is a good reference 
basis. 

b) The collision or coalescence scenario. 

Protostars are moving around in a young cluster and collisions of intermediate 
mass protostars may lead to the formation of massive stars. This interesting 
possibility and its consequences have been extensively studied by Bonnell, Zin- 
necker and colleagues (c/. Bonnell et al. 1998). 

Often in the literature (c/. Stahler 1998), it is said that the coalescence 
scenario is necessary, on the basis of the argument that the accretion scenario 
is not possible for massive stars because their high radiation pressure on the 
dust may reverse the infall and prevent the accretion. We emphasize that the 
coalescence scenario may well be important, but not for this specific reason. The 
reasons are the following ones. 

We may consider that the accretion rate M accr is given by 

Maccr ~ C S 3 /G (1) 

as resulting from the ratio of the Jeans mass by the free fall time. For typical 
temperatures of 10-100 K, as observed in molecular clouds, a value of M accr ~ 
10~ 5 Mq • yr -1 is obtained from relation (|lj). It is true that for such values of 

M accr , the effects of the radiation pressure P ra d can reverse the infall. However, 
this argument ignores the role of turbulent pressure in the clouds. Indeed, 
high turbulent velocities have been observed in regions of massive stars (c/. 
Tatematsu et al. 1993; Caselli & Myers 1995; Nakano et al. 1995; Nakano et al. 
2000). If we account for the turbulent velocities in equation (|l]), there is more 
support in the clouds which are thus denser and one gets much higher accretion 
rates of the order of 

M accr « 1O- 3 M - yr- 1 . 

For such a high accretion rate, the momentum of the infalling material is higher 
than the momentum in the outgoing radiation of a massive star and the accretion 
can not be reversed (see below, § pMj ). Thus, the basic argument often used 
against the accretion scenario is invalid, since it does not account for the fact 
that the accretion rates are very high. 
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This does not mean that the coalescence scenario never works. But, one 
has to be careful about the question of its timescale. In particular Elmegreen 
(2000) has suggested that "there is not enough time for a protostar to move 
around in a young cluster and to coalesce with other protostars" . Certainly col- 
lisions do occur, but the real importance of this scenario has still to be ascertain. 

c) The accretion scenario. 
The accretion scenario with constant was first proposed by Beech & 

Mitalas (1994). It is however evident that with moderate accretion rate like 
M accr 10 -5 Mq • yr _1 , we would need 10 7 yr to form a 100 Mq star. At this 
age, it would already have exhausted its central hydrogen. Thus, one needs 
accretion rate growing with time or with the stellar mass (c/. Bernasconi & 
Maeder 1996; Norberg & Maeder 2000; Behrend & Maeder 2001; McKee, this 
meeting). The various properties of these models will be discussed in section El 

3. Brief summary of the observational evidence for the formation of 
massive stars by accretion 

As shown by several recent works as well as in this meeting, there are several 
observational indications in favour of the massive star formation by accretion. 

• Massive outflows. 

There is evidence from radio observations of massive outflows ( cf. Church- 
well 1998; Henning et al. 2000) with outflow rates M accr oc L 75 where L is 
an estimate of the stellar luminosity. A fraction 15% to 30% of the global 
infall is supposed to be accreted, while the rest goes into the outflows. 
This suggests that during the formation of massive stars, the accretion 
rate may grow relatively steeply with the already accreted stellar mass. 

• Velocity dispersion. 

High velocity dispersions are observed, as for example from the CS line in 
the Orion KL Nebula (Caselli & Myers 1995). This supports the idea that 
turbulence provides enough support to the cloud to allow a high enough 
density, which then leads to accretion rates of the order of 10 -3 Mq -yr _1 . 

• Luminosity. 

The luminosity of IR sources like the Orion IRC2 K-luminosity is quite 
compatible with luminosities of L oc G A ^ — of discs having accretion rates 
of the order of 10~ 2 to 10~ 3 Mq ■ yr _1 (Morino et al. 1998; Nakano et al. 
2000). 

• Spectra. 

Also the spectral distribution of the energy emitted by hot cores in hot 
molecular clouds is compatible with discs having accretion rates as given 
above (Osorio et al. 1999). 

• Direct evidence for discs. 



4 



A. Maeder, R. Behrend 



Direct evidence for discs has been provided. An example is IRAS 20126 
+4004 from NH3 line observed by the VLA (Zhang et al. 1998; Cesaroni 
2000). The disc looks perpendicular to the molecular outflow, and a veloc- 
ity gradient is observed in the disc. Further compelling evidence of discs 
has been presented at this meeting by Zhang and by Conti. 

Interestingly enough, a prediction of the coalescence scenario is that discs should 
be destroyed by the collision, thus the visibility of discs may be an evidence 
against the coalescence. 



4. The accretion scenario: semi-empirical approach 



4.1. Generalities 

We are using various approaches, here. The semi-empirical approach which 
gives rough overall constraints, the analytical one which emphasizes the relation 



between basic parameters, and the accurate numerical one (cf. § 5.2.). 

There is a fundamental difference between low and high mass stars, from 
the point of view of star formation if we consider the free fall time 

ts = vim (2) 

and the Kelvin-Helmoltz timescale 

tKH = ITT ' (3) 

where p\ n is the average initial density of the cloud supposed to be at Jeans limit 
(i.e. i?j oc Mj). For M < 8Mq (in the numerical model of Figure [l], this limit 
turns to be around 15 M Q ), one has 



ts < *KH- 



This means that the accretion process is completed long before the central con- 
traction has initiated nuclear burning. Thus, in this range of masses, the evo- 
lution is not so far from the evolution at constant mass. For M ^ 8Mq, one 
has 

ts > *KH- 

In this case, the accretion is not yet completed when the central core has already 
finished its contraction and started nuclear burning. Thus, massive stars start 
their H-burning largely hidden in their molecular clouds, a fact consistent with 
the observations by Wood &; Churchwell (1989). 

Figure || illustrates several properties of the evolution with accretion. We 
start the evolution from a 0.7 Mq star which is accreting at the indicated rate 

M accr . The birthline is the path in the HR diagram followed by a star accreting 
at a specified rate, significant for dominating the evolution. If at some stage 
on the birthline the accretion is stopped, then the star sets on an individual, 
more horizontal, track corresponding to its mass and moving to the ZAMS. 
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Figure 1. Evolutionary tracks with accretion rate given by equation 
(6) with / = 0.5. The long and short dashed line is the ZAMS. Firstly 
the stars are fully convective, then radiative cores are developed and 
grow in size until a temporary fully radiative state is attained, the stars 
have finally a convective core. The age and mass of the stars when they 
quit the birthline are indicated in units of 10 3 yr and Mq at the end 
of the individual tracks, from Behrend & Maeder (2001) 



Certainly, in reality, accretion does not stop abruptly, but progressively with 
some exponential decline. However, if the decrease is fast enough, this will 
make no difference. The star finally reaches the ZAMS following a track not too 
different from the track at constant mass. However, the timescale on the tracks 
can be different from those with constant mass evolution. The accretion rate 
used in Figure |] is basically one half of the outflow rate given by Churchwell 
(1998) and Henning et al. (2000). 
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We notice in Figure |l] that the time required to form massive stars of dif- 
ferent masses are not very different, as for example for the models of 12 and 
15.7 Mq. This directly results from the accretion law used to construct the 
birthline. The birthline joins the ZAMS near 8Mq for low accretion rate of the 
order of 1CP 5 Mq • yr -1 (Palla & Stahler 1993), while for the higher accretion 
rate used in Figure |], it joins the ZAMS near 15 Mq . Above this critical point 
(which depends on M accr ), the birthline more or less coincides with the ZAMS. 
This means that the massive accreting stars are ascending the ZAMS, as they 
continue their accretion. This is a new kind of a pre-MS track: the upward track 
along the ZAMS of a rapidly accreting star. Of course, if the accretion rate does 
not keep high enough, we may have a progressive evolutionary displacement 
along post-MS tracks. 

There are several differences between the present models on the ZAMS, 
due to their particular history, and models which would be just started homo- 
geneously on the ZAMS (c/. Bernasconi & Maeder 1996). 1) A newly formed 
massive star with M > 40 Mq at the time it emerges from its cloud has already 
burnt several percents of its hydrogen. 2) A proper ZAMS does thus not exists 
for massive stars. 3) The part of the MS lifetime, during which the star is visible, 
is thus reduced. 4) The size of the convective cores are 10% smaller than for 
homogeneous classical ZAMS stars. 

4.2. Sensitivity of the birthline to the accretion rate 

Figure |2] shows the relation between the radius and the mass for stars on various 
theoretical birthlines. The noticeable point is that the higher the accretion rate, 
the more deuterium is available for burning, the larger the radius and the higher 
the birthline in the HR diagram. 

M accr f ==^ More D-burning Bigger R ==>■ More luminous birthline 

This sensitivity of the birthline to M accr offers an interesting constraint on the 
accretion rates. This was discussed with some details by Norberg & Maeder 
(2000). They assume accretion rates of the form 

M accr = A(M/MqY (4) 

and showed that values A ~ 10~ 5 Mq ■ yr -1 and tp ~ 1.0 — 1.5 are leading to 
the best fit of the birthline on a large set of observational data for T-Tauri stars 
and Ae/Be Herbig stars in the HR diagram. 

The lifetimes tpMS of the pre-MS phase also depend very much on the 
accretion rates. Figure [| shows the duration ipMS as a function of the exponent 
cp in expression (g). Note that in Tables 1, 2 & 3 of Norberg & Maeder (2000) 
there is misprint of the decimal point, the ages should be divided by 10. The 
rest of numbers in the tables are correct. Figure || shows that in order to have 
pre-MS lifetime of the order of 0.5 • 10 6 yr or less, it is necessary to have (p p> 1.0 
If this requirement is relaxed to 10 6 yr, one only needs </? ;> 0.6 These results 
support the idea that we need accreting rates growing relatively fast with the 
actual stellar mass. 
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Figure 2. Mass vs. radius relation for stars on the birthline. The 
curves labeled / = 0.5, 0.4, 0.3, 0.176 correspond to accretion rate 
which are equal to / time the outflow rates given by Churchwell (1998) 
and Henning et al. (2000), see equation (6). The track N&M is from 
Norberg & Maeder (2000). Points 1: central ignition of deuterium, 
3: start of shell D-burning, 5: end of central contraction; 6: central 
D-exhaustion, 7: H-ignition, 8: ZAMS 



4.3. Relation with the outflows observed by Churchwell and Henning 

Massive outflows have been identified from radio and IR observations by Church- 
well (1998) and Henning at al. (2000). The remarkable fact is that there is a 
relation between the stellar luminosity and the outflow rate M out of the form 
M out oc L - 75 . According to a paper presented by Beuther et al. at this meeting, 
the relation observed by Churchwell and Henning et al. could rather be an up- 
per envelope of the M out distribution as a function of luminosity. If we adopt a 
typical mass luminosity relation for stars on the ZAMS in the range 2 — 85 Mq , 
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Figure 3. Time needed by a star to grow from 1 to 60 Mq at a rate 
specified by expression (|j), as a function of (p with A = 10~ 5 Mq ■ yr _1 



we get 

M out « 1.5 10~ 5 (M/M ) L5 Mq • yr- 1 . (5) 

If we assume that a fraction / of the infalling material is accreted by the star 
and (1 — /) is going into the outflow, we get that a fraction / = //(l — /) of the 
outflow mass is effectively accreted, i.e. 



f M, 



out • 



(6) 



This implies that the accretion rates grow swiftly with the stellar mass, like 

Figure | 



4.2 



M . This result is quite consistent with what we have found in 
shows numerical models by Behrend & Maeder (2001) made with equation (6) 
and various values of /, where M out is taken from data of Churchwell (1998) and 
Henning et al. (2000). We see that models between those with / = 0.5 and 0.176 
beautifully reproduce the upper envelope of the observations. A value / = 0.5 
means that one third of the infalling material is effectively accreted, / = 0.176 
means 15%. These values are those suggested respectively by the models by Shu 
et al. (1998) and by the observations by Churchwell (1998). 

Table |l] shows some properties of the models by Behrend & Maeder (2001). 

The disc luminosity is taken as Ldisc = G M ^ llflow with Minnow = M out + M accr 
(c/. Hartmann 1998). Of course there is a proportionality factor which is un- 
certain and thus these values for L^[ sc are only indicative. Table |l] shows that 
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Figure 4. Birthlines in the HR diagram for various values of / from 
Behrend & Maeder (2001). Values / = 0.5 and 0.176 correspond to 
1/3 and 15% of the infalling material accreted. A minimum value of 
10~ 5 Mq -yr -1 is taken for low mass stars. The observations are those 
collected by Norberg & Maeder (2000); the birthline proposed by these 
authors is indicated by N& M. 

the last stages of massive star accretion go very fast. For example, in Table |] 
the time for the evolution from 15 to 80 Mq is 4.7 10 4 yr. For / = 0.5, this time 
would be only 1.5 10 4 yr. Of course, for different values of <p as in equation (|j), 
we would have different timescales as suggested by Figure |j| As we do not 
know exactly what is the value of (p and /, this means that the lifetime for the 
fast phases of massive star accretion, i.e. from 15 Mq to 80 Mq are still rather 
uncertain, somewhere between a few 10 4 yr and a few 10 5 yr. 
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Table 1. Properties of stars on the birthline with / = 0.176: age in 
units of 10 3 yr, luminosity, mass, effective temperature, radius, relative 
upper radius of the convective zone, ratio of disc to star luminosities. 
From models by Behrend & Maeder (2001). 
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5. Models and discussion 

5.1. Physics feasibility of the high accretion rate 

The accretion rates given by equations (|j) and (^) are remarkably well located 
in the stable zone derived by Wolfire & Cassinelli (1987). For rates lower than 
a certain limit, the momentum in the radiation pressure is sufficient to prevent 
the accretion. For rates higher than some limit, the luminosity created by the 
shock is supra-Eddington. It has been shown by Nakano (1989) that for a 
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non-spherical collapse the stability conditions are much less severe than those 
given by Wolfire h Cassinelli. In particular, Nakano (1989) suggests that, even 
for the most massive stars, radiation cannot stop the accretion process and as 
a consequence he suggests that the maximum stellar mass is determined by 
fragmentation rather by the accretion process. Certainly, both fragmentation 
and accretion are essential in shaping the IMF and determining the maximum 
mass. Indeed, if most of the infalling material is diverted in the outflows, this 
implies that only a tiny fraction of the initial fragmentats finally contributes to 
the IMF. 

We also point out that the various accretion models do not account for the 
anisotropy of the radiation field of massive rotating stars, which would favour 
anisotropic accretion in pre-MS stages and anisotropic mass loss in post-MS 
phases (cf. Maeder & Desjacques 2001). Indeed, the lower T e g in the equatorial 
regions of a rotating massive pre-MS object should very much favour the accre- 
tion on these objects. Joined to the above results of Nakano (1989), this supports 
the idea that accretion at heavy rates of the order of M accr ~ 10~ 3 Mq • yr _1 is 
quite possible theoretically. 

5.2. Dependence of the accretion rate on the central mass 

The location of the birthline, the timescales and the observations of the outflows 
support accretion rates M accr ~ 10 -5 (M/Mq) 1 ^ with (p ^ 1 as shown above. 
Globally, we may understand a relation of this kind since for larger masses the 
central potential is much larger and attracts more matter. The accretion must 
go like 

Maccr OC ^ OC M (Gp)V2 a (M/i?) 3 / 2 
Iff 

where M is the mass of the initial cloud of average density p, tg is the free fall 
time. The radius of the initial subcondensation may vary with their mass in 
different ways according to different authors. Larson (1981) has suggested for 
small condensations (Larson's scaling) 

'"it 

If so, R oc M 1 / 2 and thus M accr oc M 3 / 4 . In this meeting, Johnstone (2001) 
has suggested that the density of the subcondensations is constant (Johnstone's 
scaling) 

p oc const. 
In this case, R oc M 1 / 3 and we have 

M accr oc M. (7) 

Interestingly enough, Johnstone's scaling implies accretion rate growing almost 
linearly with the mass as suggested by the observations quoted above. The above 
arguments are rather schematic and ignore many effects like rotation, accretion 
disc, outflows, etc. Nevertheless, the general trend they show must be roughly 
correct. Very detailed numerical models are now in progress in Behrend's thesis. 
With finite elements grids, which allow us to follow more closely the interesting 
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parts of the models, the accretion with rotation is modelized, the inflow of the 
matter through the disc followed as well as the central accretion and outflows. 
These models predict the accretion rates and should be combined with interior 
models of the central body. 

5.3. Formation of binaries 

An apparent difficulty of the accretion scenarios (it is sometimes presented as 
objection) is the question on how massive binary stars do form. The question is 
especially important because the fraction of binaries and multiple stars is very 
high as shown by Zinnecker (this meeting). Here, we would like to suggest some 
possibilities of preferential binary formation among the OB-stars. 

Let us consider a molecular cloud, which is collapsing to give birth to a 
young star cluster. In the process of cloud fragmentation, some cores may be 
bound gravitationally, some do not. Each core, bound or not, will accrete matter 
from its own appropriate Jeans radius. However, the bound or multiple cores will 
also, in addition to their current accretion, experience a Bondi-Hoyle accretion 
particularly if their Bondi-Hoyle radius b 

GM 

b = 2— (8) 

v 

is greater than the semi-major axis of the supposed double system. In this 
case, the velocity v appearing in equation (^|) is some combination of the orbital 
velocity in the double system and of the velocity of the double system in the 
cluster. Thus, double and multiple systems have a higher potential reservoir of 
matter for accretion and their accretion rates M accr should also be larger than 
that for single stars. 

This means that double and multiple cores, by their larger reservoir and 
higher M accr will have more chance to move to the top of the mass scale (the 
IMF) and they will do it faster than the other single objects. This may thus 
lead to a higher number of double and multiple systems among OB-stars in 
young clusters. Since mass is accreted on the protobinary system, the semi- 
major axis will decrease and the orbital periods will also decrease. It is to be 
examined whether an orbital separation corresponding to the initial tidal radius 
within the cluster finally leads to periods P < 10 d, as so frequently observed in 
massive O-stars. 

In conclusion, we see that many observations and theoretical considerations 
effectively support the accretion scenario for massive stars. The high frequency 
of binaries among O-type stars might even be an additional argument in favour 
of the accretion scenario. 
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